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We systematically study gapless topological phases of (semi-) metals and nodal superconductors 
described by Bloch and Bogoliubov-de Gennes Hamiltonians. Using K-theory, a classification of 
topologically stable Fermi surfaces in (semi-) metals and nodal lines in superconductors is derived. 
We discuss a generalized bulk-boundary correspondence that relates the topological features of the 
Fermi surfaces and superconducting nodal lines to the presence of protected zero-energy states at the 
boundary of the system. Depending on the case, the boundary states are either linearly dispersing 
(i.e., Dirac or Majorana states) or are dispersionless, forming two-dimensional surface flat bands or 
one-dimensional arc surface states. We study examples of gapless topological phases in symmetry 
class AIII and Dili, focusing in particular on nodal superconductors, such as non-centrosymmetric 
superconductors. For some cases we explicitly compute the surface spectrum and examine the 
signatures of the topological boundary states in the surface density of states. 
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I. INTRODUCTION 

The recent discovery of topological electronic phases in 
insulating materials with strong spin-orbit coupling [T]-[6] 
has given new impetus to the investigation of topological 
phases of matter. Topological materials, such as the inte- 
ger quantum Hall state and the spin-orbit induced topo- 
logical insulators, are characterized by a nontrivial band 
topology, which gives rise to protected exotic edge (or 
surface) states. Many interesting phenomena, including 
magneto-electric effects [7] and the emergence of local- 
ized Majorana states, [8] have been predicted to occur in 
these systems. These phenomena could potentially lead 
to a variety of new technical applications, including novel 
devices for spintronics and quantum computation. 

Besides the topological insulators and the integer quan- 
tum Hall state, which have a full bulk gap, there are also 
gapless phases that belong to the broad class of topolog- 
ical materials, such as, e.g., (semi-)metals with topologi- 
cally protected Fermi points and nodal superconductors 
with topologically stable nodal lines. These gapless topo- 
logical phases also exhibit exotic zero-energy edge (or 
surface) states with many interesting properties. These 
boundary states may be linearly dispersing (i.e., of Dirac 
or Majorana type), or dispersionless, in which case they 
form either two-dimensional surface flat bands or one- 
dimensional arc surface states. Notable examples of gap- 
less topological materials include, among others, [S'-TT 
graphene, [12-14 (ia,2_^2-wave superconductors, [15-,.-,! 7^ 
the A phase of superfluid ^He, [l¥l |19] and nodal non- 
centrosymmetric superconductors. [2QH28] 

The topologically stable Fermi points and supercon- 
ducting nodal structures in the aforementioned materi- 
als can be viewed, in a sense, as momentum-space de- 
fects, that is, as momentum- space analogues of real- 
space topological defects. In other words, the nodal 
points in dx2_y2-\^me superconductors, the Fermi points 



in graphene, and the nodal points in ^ He A can be inter- 
preted as momentum-space point defects, i.e., as vortices 
and hedgehogs, respectively. The nodal lines in non- 
centrosymmetric superconductors, on the other hand, 
correspond to momentum- space line defects, i.e., vortex 
lines. Similar to real-space defects, the stability of these 
Fermi points, nodal points, and nodal lines is guaranteed 
by the conservation of some topological invariant, i.e., 
e.g., a Chern or winding number. 

In this paper, building on previous works, |T8l [27]- 
l38l [40] we derive a classification of topologically sta- 
ble Fermi surfaces in (semi-) metals and nodal lines in 
superconductors using K-theory arguments. [43] As it 
turns out, the presence of discrete symmetries, such as 
time-reversal symmetry or particle-hole symmetry, plays 
a crucial role for the classification of gapless topological 
phases, a fact that has not been emphasized previously. 
The appearance of protected zero-energy states at the 
boundary of gapless topological phases is discussed and 
it is shown that the existence of these boundary states 
is directly linked to the topological stability of the Fermi 
surfaces (superconducting nodal lines) in the bulk via a 
generalized bulk-boundary correspondence. In particu- 
lar, we demonstrate that gapless topological phases sat- 
isfying the so-called chiral symmetry necessarily support 
zero-energy surface flat bands. 



II. LOCAL STABILITY OF FERMI SURFACES 

The classification of topologically stable Fermi surfaces 
in terms of K-theory is closely related to the classification 
of topologically stable zero modes localized on real space 
defects. 



In Sec. HA 



we will therefore first review the 
stability of localized gapless modes on topological defects, 
before discussing the classification of topologically stable 



Fermi surfaces in Sec. II B To denote the dimensionality 



2 



of the Brillouin zone (BZ), the Fermi surfaces, and the 
real space defects we use the notation: 

dBZ = (total spatial dimension) 

= (total dimension of the BZ), 

dk = (codimension of a Fermi surface) — 1 

_ / # of parameters characterizing a surface \ 
I surrounding a Fermi surface in the BZ i ' 

dr = (codimension of a real space defect) — 1 

_ / # of parameters characterizing a surface \ 
~ y surrounding a real space defect J ' 

In other words, the dimension of a Fermi surface and 
a real space defect are c^bz — dk — 1 and c^bz — dr — 1^ 
respectively. 

A. Real space defects 

In this subsection, we review the stability of local- 
ized gapless modes on topological defects from the K- 
theory point of view. [37| [44] To that end, let us con- 
sider the topology associated with gapped Hamiltonians 
?^(r, /c), where k are the bulk momenta and r are real 
space parameters characterizing the defect. That is, r are 
the coordinates parametrizing the surface that encloses 
the defect in question. For instance, a line defect in a 
three dimensional system is described by the Hamiltonian 
n{r,k) = kiji + /C272 + ksjs + mi{x,y)j4 + m2{x,y)j3, 
where mi{x,y) = xj^x^ + y'^ , m2{x,y) = yj \/x^ + y^ , 
and 7's are 4x4 anti-commuting matrices. In this case, 
^1,2,3 ^ k and mi ^2 ^ ^- The stable equivalent classes of 
Hamiltonians H(r, k) can be described by the K-group 

K^{s\dx,d2), (1) 

where s represents one of the Altland-Zirnbauer symme- 
try classes given in Table [!) and F = C (R) for the com- 
plex (real) Altland-Zirnbauer symmetry classes. d\ and 
d^ represent the dimension of r and /c, respectively. 

For condensed matter systems defined on a lattice, the 
BZ is a (i2-dimensional torus, k G T^'^ ^ and r G 5*^^ where 
5'^^ is a -dimensional sphere surrounding the defect in 
real space. If we are interested in "strong" but not in 
"weak" topological insulators and superconductors, we 
can take k G 5''^^. Furthermore, it turns out it is enough 
to consider (r,k) G 

An underlying strategy in Ref. 37 is that (i) we move 
apart from the defect (far enough), (ii) adiabatically go 
around the defect, and (iii) monitor the wavefunctions of 
the family of Hamiltonian H(r, /c). From this Hamilto- 
nian, one can define a K-theory charge for the system. 
One of the main results in Ref. |37l is the connections of 
K-groups among different systems 

K^{s\dx,d2 + 1) = K^{s - \\dx,d2), 

K^{s;di + 1,^2) = K^{s^l;di,d2). (2) 
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TABLE I: K-group Kw{s;d2,di) from Ref. [37]. 



In particular, by combining these two relations, 

K^{s; di + 1, ^2 + 1) = K^{s; 6/1, 6^2), (3) 

which shows physics depends only on 

S = d2 - di. (4) 

These relations ([3| can be obtained by explicitly 
connecting two Hamiltonians which belong to different 
symmetry classes and different position or momentum 
dimensions. Since the connection is smooth, the two 
Hamiltonians are topologically equivalent, i.e., they are 
in equivalent classes. 

One of the key issues in their analysis is to know 
how each anti-commuting matrix and its coefficient in 
a non-interacting Hamiltonian transforms under the dis- 
crete symmetries. For instance, suppose we decompose 
a Hamiltonian H(r, k) with position type Dirac matrices 
7^ and momentum type Dirac matrices 7^, 

n{T,k) = R^{r,k)-i^^ K\r,k)-i,. (5) 

Here, {7/i,7iv} = ^^^iv, {lulj} = 2^ij, {7/i, 7i} = 0, and 
for time reversal symmtery T, we require 

[7p,T] = {7i,r} = (6) 

while for particle-hole symmetry P, 

{l^,P} = h,P]=0. (7) 

R and K should transform in the same way as position 
and momentum. If a Hamiltonian contains 6 + 1 position 
like matrices 7^ and a momentum like matrices 7^, then 
the Hamiltonian belongs to a symmetry class 

a — b = s mod 8. (8) 
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TABLE II: K-group Kw{s; di, d2) where di and c/2 are rel- 
evant dimensions in real and momentum space, and classifi- 
cation of (i) bulk insulators/superconductor, (ii) bulk insu- 
lators/superconductor with defect, (iii) Fermi surfaces, (iv) 
Fermi surfaces with defect. 



^From the K-group Ky{s; di^d2) and its shifting prop- 
erties (Table [l|), one can reproduce the periodic table of 
topological insulators and superconductors by taking 



di =0, d2 = d] 



BZ, 



d] 



BZ- 



(9) 



Teo and Kane in Ref. [37 use K^{s]di^d2) to predict 
and classify the appearance of zero-energy modes in the 
presence of a topological defect in real space. Here, we 
take 



di = dr 



do 



^BZ, 



5 = d^ 



BZ 



(10) 



When the corresponding K-group K^{s]di^d2) is non- 
trivial, there is a zero energy mode(s) localized on the 
topological defect. 



B. Fermi surfaces (Momentum space defects) 

The analysis of Ref. '37', which we have reviewed above, 
can be extended to study the stability of Fermi surfaces. 
In the presence of a q{= ^bz — dk — l)-dimensional Fermi 
surface in the BZ, we move apart from the Fermi sur- 
face and monitor the wavefunctions of the family of the 
Hamiltonian 1-L{k) in momentum space where k G S^^ = 
gd^z-q-i ^ gdBz = BZ is the hypersurface in the BZ 
surrounding the Fermi surface. It then turns out we can 
take 



di =0, d2 = dk, S = dk. 



(11) 



I.e., we look at Ky{s;5 = dk). The proposed classifica- 
tion (or "periodic table") of topologically stable Fermi 
surfaces is presented in Table [l[ 

To illustrate the discussion above, we list examples of 
stable (and also unstable) Fermi surfaces. We consider 
a single particle Hamiltonian in momentum space H(/c), 
where k = {ki^k2, ■ ■ ■ ,kd) denotes the d-dimensional mo- 
mentum coordinate. The Fermi surface is defined by 
1-L{k) = 0. More realistic examples will be given in the 
later sections. First, let us consider class A. This is the 
symmetry class with no constraints on the Hamiltonian. 
We here list the examples of class A Hamiltonians which 



have a Fermi surface, together with the dimension of the 
Fermi surface (2nd column): 



Hamiltonian 

n{k) = ki 

1-L{k) = kiCFi 
1-L(k) = kiCFi 
1-L{k) — kiai 



k2(72 

k2(J2 - 
- k2(y-2 



ksCTs 
- ksas 



dimension(g') 
d-1 
d-2 

d-3 (12) 
d- 4 



(Here, cri,2,3 are the Pauli matrices, 0^1,2,3 and /3 are the 
Dirac matrices.) The example with q = ^bz — 1 de- 
fines the Fermi surface /ci = 0, which is clearly stable 
against any perturbation to the Hamiltonian, as Horava 
predicted. [30] Similarly, the example q = c^bz — 3 de- 
fines the Fermi surface /ci = /c2 = /C3 = 0, which is also 
stable against any perturbation. On the other hand, ex- 
amples with (^dk -\- I = dBZ — q) even are not stable. In 
class A, there is no spectral symmetry (i.e., the chem- 
ical potential is not pinned at /i = 0). Therefore, for 
example, what appears to be a stable Fermi "point" for 
1-L{k) = kiCTi + k2(T2 + k^^cj^, when d^z = 3 can actually 
represent a stable Fermi surface if the chemical potential 
is included. 

Next, let us consider class AIII. Remember that a 
Hamiltonian in symmetry class AIII satisfies, by defini- 
tion, one symmetry condition, 



{n{k),Us} = o, 



(13) 



where Us is an arbitrary unitary matrix. We here list the 
examples of class AIII Hamiltonians which has a Fermi 
surface, together with the dimension of the Fermi surface 
(2nd column): 



Hamiltonian 
1-L{k) — kidi 
1-L{k) = kiCFi - 
1-L(k) = kiai 
1-L(k) = kiai 



k2(72 

- k20L2 

- k20i2 



ksas 
ksas 



dimension (g') 
d-1 
d-2 

d-3 (14) 

k^(3 d-4 



^From these general K-theory arguments, we can ob- 
tain a classification of topologically stable single gapless 
mode (Table [m). 

Note that the above analysis can also be applied to 
Hamiltonians living in an extended parameter space, i.e., 
Hamiltonians that are parametrized by the momentum 
coordinates and some external control parameters, such 
as e.g. mass terms m^. The topological arguments can 
then be used to predict the existence of finite regions of 
gapless phases in the topological phase diagram. 



C. Comments 

It should be stressed that the above stability criterion 
is for a single Fermi surface; so far we are primarily in- 
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TABLE III: (Single Fermi surface) Classification of single 
Fermi surface in metals and superconductors as a dimension 
dBz — c/fc — 1 of the Fermi surfaces and symmetry class of the 
single Fermi surface in c/bz = 2,3 spatial dimensions. 



terested in S^'^ which is surrounding one Fermi surface. 
In many situations, there can be multiple Fermi surfaces 
that are located in different regions in the BZ. While 
in such cases S^^ can be extended such that it encloses 
all the Fermi surfaces of interest to study their stabil- 
ity, it is important to notice that there is an intertwining 
role played by lattice effects and discrete symmetries, on 
which we make comments below. 

a. Fermion doubling It should be noted that in lat- 
tice systems, there may be the fermion doubling problem; 
In lattice systems, Fermi surfaces may be accompanied 
by their partner once one imposes discrete symmetry con- 
ditions. I.e., Fermi surfaces are always created pairwise. 
This is to some extent expected since for any lattice sys- 
tems, there should be a way to gap Fermi surfaces, by 
nesting, say. If this is the case, if S^^ is extended so 
that it encloses all the Fermi surfaces of interest, the K- 
theory charge should always be zero. In other words, 
the above consideration for the local stability is useful, 
when we confine ourselves to situations where we switch 
off "nesting" of all kinds. 

b. Effective symmetry classes As we consider each 
Fermi surface in the BZ separately, one should also be 
aware that the determination of the relevant symmetry 
class needs a careful treatment. For instance, suppose 
we consider a time reversal symmetric system which has 
two Fermi points related to each other by time rever- 
sal. The full system is time reversal symmetric, and we 
can consider in principle a topological number associated 




FIG. 1: (Color online) Illustration of integration path defor- 
mation. (/ci|| , /c2|| ) are the momenta in the boundary direc- 
tions and ki± is the momentum in the orthogonal direction 
to the boundary. 



with this symmetry class to study the stability of those 
gapless modes. On the other hand, each of the Fermi 
points itself is not time reversal symmetric. So the topo- 
logical number associated with this single Fermi point is 
not the one defined for time reversal symmetric Fermi 
points, and the single Fermi point belongs to a differ- 
ent symmetry class. So even though we are looking at 
the same system, we need to consider different symme- 
try classes depending on whether we are interested in the 
single Fermi point or the total system. 

c. Alternative criterion from higher dimensional 
topological insulators and superconductors Instead of 
considering the local stability of the Fermi surfaces as 
above, one could construct topologically stable Fermi 
surfaces from the boundary modes of topological insu- 
lators and superconductors as follows. The classification 
of topological insulators and superconductors includes, 
implicitly, a classification of stable gapless states; these 
are the states that appear at the boundary of the topo- 
logical bulk, and the existence of stable gapless states at 
the boundary and the topological bulk has a one-to-one 
correspondence. The stable gapless boundary modes of 
(i-dimensional topological bulk is a stable Fermi surface 
in d— 1 dimensions. In fact, the bulk topological invariant 
is directly related to the topological charge defined for the 
boundary Dirac fermion modes, as discussed in Ref. 1391 
Furthermore, we argue that such gapless Fermi surface 
can be embed in the d' dimensional BZ, where d' can dif- 
fer from d {d' > d). Even when d' 7^ the Fermi surface 
remains stable. This suggests that the classification of 
locally stable Fermi surface can be obtained by a simple 
dimensional "shift" of the periodic table of topological 
insulators/superconductors. The result is summarized in 
Table [El 

In this case, the symmetry of the gapless state means 
that of the total system, rather than a single Fermi point. 
Note also that we should again note the Fermion dou- 
bling problem. The gapless boundary modes of ^bz- 
dimensional topological bulk is an absolutely stable Fermi 
surface md—1 dimensional BZ; they never go away, as far 
as the bulk topological character is not destroyed. How- 
ever, when we embed them in the d^^r^ / dez dimensional 
BZ, such Fermi surfaces must be accompanied by their 
partner. They can be pair annihilated. 
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TABLE IV: (Total system) Classification of Fermi surfaces 
in metals and superconductors as a dimension c/bz — c/fc — 1 of 
the Fermi surfaces and symmetry class of the total momentum 
space in c^bz = 2,3 spatial dimensions. 



D. Topological invariants and surface flat bands 

As we mentioned, gapless topological insulators and 
superconductors can exhibit surface flat bands. In this 
section, we systematically classify stable flat bands from 
K-theory point of view. Let us first consider a Fermi 
ring in a three-dimensional BZ. In this case, the hyper- 
surface surrounding the Fermi ring is one dimensional: a 
circle for a lattice model or a line for a continuum model. 
If the topological number associated with the induced 
Hamiltonian on the line is non-trivial, the Fermi ring 
is stable. Now we deform the circle so that it pierces 
the Fermi ring as shown in FIG. [l] This deformation 
is analogous to that of a contour deformation in com- 
plex analysis. Let us denote the momenta parallel to the 
boundary as {ki\\,k2\\) and orthogonal to it as ki±. The 
final circle extends along kij_ for a fixed (/ciy , A:2||). If a 
circle/line for a given pierces only a Fermi ring 

with non-trivial K-theory charge, there must be a zero 
mode at this momentum value on the boundary in real 
space. If a circle/line for a given (^f||,^|||) pierces only 
Fermi rings with both positive and negative K-theory 
charges and the net charge is zero, or does not pierces 
a Fermi ring at all (contractible), then there will be no 
topologically protected gapless modes at that value of 
{k^^k^). Since this argument applies for each value of 

(^f||' ^2||)' ^here will be a two-dimensional region of gap- 
less modes (i.e., a flat band) in the surface state. Next, 
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TABLE V: The dimensions of the Fermi surfaces and the flat 
bands 



we consider a Fermi ring in a four-dimensional BZ. In 
this case, the topological charge associated with it is de- 
fined on a two-dimensional sphere or torus surrounding 
it. We parametrize it by {ki±^k2±) for a fixed (/ciy , /c2||). 
Then, again, the gapless mode appears for each value of 
(^i||'^2||)' which leads to a two-dimensional flat band. 
With a similar argument, we conclude that a Fermi ring 
always gives a two-dimensional fiat surface band. 

We also look at a Fermi point in a three-dimensional 
BZ. A Fermi surface is surrounded by a two-dimensional 
sphere and this surface is deformed to a two-dimensional 
torus. We parametrize it by (/cix,A^2±) for a fixed ki\\. 
We assign a topological charge for each value of /ciy. So 
the dimension of the surface flat band is one. For a Fermi 
point in a four-dimensional BZ, the surrounding surface 
is three dimensional and replace the parametrization to 
(^i±,^2±,^3±) for a fixed /ciy. Again, we obtain a one- 
dimensional flat band parametrized by ki\\. The result is 
summarized in Table fV} By these observations, it is now 
clear that the dimension of surface flat band is given by 
the dimension of the Fermi surface +1. 

Before closing this section, let us comment on the 
case dk = 0. This includes a Fermi surface in a three- 
dimensional BZ. In this case, the meaning of a hyper sur- 
face surrounding the Fermi surface is not clear: a Fermi 
surface divides the total BZ into two regions and there 
is no way to surround the Fermi surface only inside the 
momentum space. One may be able to extend the space 
to a frequency-momentum space and surround the Fermi 
surface in the frequency direction. [29 Or one may stay 
in momentum space and consider a Hamiltonian induced 
on points in both sides of the Fermi surface. However, in 
either case, there is no way to project the Hamiltonian or 
the Fermi surface onto a boundary. In this sense, there 
is no gapless mode or zero-energy flat band when dk = 0. 



III. PROTECTED SURFACE STATES IN 
NODAL TOPOLOGICAL SUPERCONDUCTORS 

In order to make the argument of the previous sec- 
tion more concrete, we consider the case of nodal lines 
in three-dimensional topological superconductors with 
time-reversal symmetry. To characterize the topological 
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properties of the BdG Hamiltonians, we introduce three 
different topological invariants. We shall keep the analy- 
sis as general as possible, such that it may be applied to 
arbitrary superconducting systems. 



A. Topological Invariants 

Let us consider a general time-reversal invariant BdG 
Hamiltonian H{k). This could be either a supercon- 
ductor belonging to class Dili, AIII or CI. Remember 
that time-reversal invariant BdG Hamiltonians posses a 
chiral symmetry, which is realized as a combination of 
time-reversal and particle-hole symmetry. I.e., there is a 
unitary matrix S which anticommutes with the Hamil- 
tonian. To characterize the topological properties of 
H{k) we introduce three different topological invariants. 
ISEllaSlllSHiH]. The topological properties of H{k) (or 
the quantum ground state) is encoded in the spectral 
projector Q{k)^ 



Q{k) = J2\Mk)){ua{k)\, 



(15) 



where \ua{k)) is the eigenstate of H{k)^ and the summa- 
tion extends over the filled states. (To be more precise, 
P{k) = [1 — Q{k)]/2 is the spectral projector.) The set 
of filled states and associated wavefunctions over the BZ 
form a vector bundle. When the bundle is "twisted" , in 
a sense that is characterized by a non- vanishing topolog- 
ical invariant (Chern invariant, say), the ground state is 
in a topologically non-trivial phase. 

Since time-reversal invariant BdG Hamiltonians pos- 
sess a chiral symmetry, there is a natural grading in which 
the eigenequation H{k)'^^{k) = ±Aa(/c)^J(^) ^^kes off- 
diagonal form 



D{k) 



(16) 



where a = 1, . . . , 2N is the combined band and spin in- 
dex. We assume that there is a spectral gap around zero 
energy (at least in some part of the BZ) with \Xa{k)\ > 0, 
and for definitiveness we choose Xa{k) > for all a. In 
analogy to Eq. (15) in this case, the natural object to 



look at is the off-diagonal component of Q{k) 
qik) = J2Mk)){vaik)\, 



(17) 



where the eigenfunctions (x^)'?^) can be obtained from 
the eigenvectors of DD^ or D^D 



A 



(18) 



The eigenvectors Va are taken to be normalized to 

1, for all a (here, the index a 



one, I.e., ulu 



is not summed over). The g'-matrix is the off-diagonal 
block of the spectral projector (see below). 



Observe that The eigenstates of D^D follow from the 
eigenstates of DD^ via 



(19) 



Va =MaD^Ua, 

with the normalization factor Ma- Using Eq. ([l8| one 
can check that Va is indeed an eigenvector of D^D, 

D^DVa = D^D{MaD^Ua) = K^ID^u^ = A^w,, (20) 

for all a. The normalization factor Afa is given by 

uiDD^u, = Xluiua = Xl J^, = ^, (21) 
for all a. It follows that the eigenfunctions of H{k) are 



(22) 



With this, the projector P{k) onto the filled Bloch states 
becomes 



-^i:(:;)(»i-"i) 



(23) 



hN 

2\0 hNj 2 

Finally, we obtain for the flat band Hamiltonian Q, which 
is defined as Q = I4JV — 2P , 



Q{k) 



q{k) 
q\k) 



q{k) = Y,u,{k)vl{k)D{k) = Y^Ua{k)u\{k)^^. 



(24) 



1. Topological winding number 



The integer- valued topological invariant characterizing 
topological superconductors is now simply given by the 
winding number of q{k). It can be defined in any odd 
spatial dimension. In three dimensions we have 

= /bz 24^'"^'^ [iq-'d,q){q-'d.q){q-'d,q)] , 

(25) 



and in one spatial dimension it reads 



(26) 



Alternatively, it is also possible to define the winding 
number in terms of the unflattened off-diagonal block 



7 



D{k) of the Hamiltonian. For example, for the winding 
number in one spatial dimension this reads 

1^1 = ^^J^JkTT[D-^dkD-{D^}-^dkD^ 

= ^Im [ dkTi [dk In D] . (27) 
27r Jbz 

In the nodal superconducting phases the winding num- 
ber u is no longer quantized. However, we can consider 
H{k) restricted to ID loops in reciprocal space and define 
a topological number in terms of a ID momentum space 
loop integral to characterize the topology of the gapless 
phases. We observe that H{k) confined to a generic mo- 
mentum space loop no longer satisfies TRS nor PHS, but 
it still obeys chiral symmetry S. Hence, H{k) restricted 
to a loop in the BZ belongs to symmetry class AHI [33 
and its topological characteristics are described by the 
ID winding number 



Nc = £dlTv[q-\k)Viq{k)], 



(28) 



where the integral is evaluated along the loop jC in the 
BZ. Observe that for any closed loop C that does not 
intersect with gapless regions in the BZ, is quantized 
to integer values. If jC is chosen such that it encircles a 
line node, then Njc determines the topological stability 
(i.e., the topological charge) of the gapless line [27| l28]. 



2. Z2 Invariant 

In this section we compute the Z2 topological invariant 
for symmetry class Dili in d = 1 and d = 2 spatial dimen- 
sions. It is most convenient to perform this derivation 
in the basis (16), in which the AN x AN Bogoliubov-de 



Gennes Hamiltonian takes the form 



H{k) = 



D{k) 
D^k) 



D{k) = -D'^i-k). (29) 



In this representation, the time-reversal symmetry oper- 
ator is given by T = K-Ut = K. i(J2 ® ^2N and the flat 
band Hamiltonian reads 

Q{k) = J^^^^ "^^^^ j , q{k) = -q^i-k). (30) 

The presence of time-reversal symmetry allows us to de- 
fine the Kane-Mele Z2 invariant [H [49H53] , 



w = i[ 



K 



Pf [w{K)] 
v/det [w{K)] ' 



(31) 



with K a time-reversal invariant momentum and Pf the 
Pfaffian of an anti-symmetric matrix. Here, w{k) denotes 
the "sewing matrix" 



Warn = (u+(-fc)|r<(fc)), 



(32) 



where a, 6 = 1, ... , 2N and u^{k) is the a-th eigenvector 
of Q{k) with eigenvalue ±1. The Pfaffian is an analog 
of the determinant that can be defined only for 2n x 2n 
anti- symmetric matrices A. It is given in terms of a sum 
over all elements of the permutation group S2n 



-l),cr(2i)- 



crGS'2. 



Due to the block off-diagonal structure of Eq. (30) a 



set of eigen Bloch functions of Q(k) can be constructed 
as |6] 



\ut{k))n = 



V2 [ ±qHk)7 



or, alternatively, as 



\utik)h 



_ J_ / ±q{k)r 
~ V2\ ria 



(33) 



(34) 



where are 2N momentum independent orthonormal 
vectors. For simplicity we choose {na)h = ^ab- In passing, 
we note that both \uf{k))]si and \uf{k))s are well-defined 
globally over the entire Brillouin zone. To compute the 
Z2 topological number we choose the basis \u^{k))^. 
Combining Eqs. (32) and (33) yields 



q\k)ni, 



(35) 



In the second last line we used Eq. (30), i.e., q{—k) = 
—q^ik). In conclusion, the Z2 topological number in spa- 
tial dimensions d = 2 and d = 1 is given by 



W 



n 

K 



Vi\f{Kj\_ 
\/det [q{K)] ' 



(36) 



where K denotes the four (two) time-reversal invariant 
momenta of the two-dimensional (one-dimensional) Bril- 
louin zone. 

In nodal superconductors, the Z2 invariant is not well- 
defined in the bulk. However, we can consider H{k) re- 
stricted to a time-reversal invariant (TRI) loop which 
is mapped onto itself under k ^ —k. In that case we 
obtain a ID Hamiltonian satisfying both TRS and PHS 
(i.e., belonging to symmetry class Dili). The topological 
properties of such a ID system are characterized by the 
following Z2 invariant 



K 



\/det [q{K)] 



(37) 
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where K denotes the two TRI momenta on the loop C. 
Note that Wc is either +1 or -1 for any TRI loop that 
does not cross gapless regions in the BZ. 



Nodal topological superconductors with 
conservation 





Here, we study a lattice version of a class AIII topo- 
logical superconductor. This model is equivalent to the 
polar state of ^He. The BdG model Hamiltonian is given 



by ^ = Efc ^i^{k)^k where = (c 



and 



n{k) 




^t,k —^k - OLk . 



(38) 



with £k = t(cos kx + cos ky + cos kz) — /J^^ ak = a [sin kz] , 
Aq [sin kz]. The eigenvalues are given by 



t,k 



and A 

Xk = \l {^k + <^/c)^ + (As + At,/c)^. Observe that the z- 
component of the spin is a conserved quantity in this 
model. The Hamiltonian anti-commutes with a2 and thus 
it belongs to class AHI. A two-layer version of this 
model might be realized in the pnictide superconductor 
SrPtAs.|54, 55 



Model (38) with At^k = Ao[sin/ca;] has one nodal ring 
which lies in the (x, ?/) -plane, centered around the T point 
of the BZ. The nodal ring is described by 

kz = 0^ kx = ^ arccos [ji/t — 1 — cos ky] . (39) 

The topological properties of the nodal lines are given by 



1^1 



where 



-ak -Sk + i {At,k + As) . 



(40) 



(41) 



For simplicity, we set A^ = in what follows. Fig. |2] 
shows the winding number for the (001) surface. Alter- 
natively the topological property of the nodal line can be 
expressed in terms of the g- matrix = '^^a^a^k^ 
where Ua are the eignevectors of D^D. 



C. Non-centrosymmetric superconductors 

a. BdG Hamiltonian In this section, we study a sim- 
ple theoretical model of time-reversal symmetric, non- 
centrosymmetric superconductors (BdG Hamiltonians) 
with only single band (per spin). The relevant symme- 
try class is DHL In addition, another important input 
for physics of non-centrosymmetric superconductors is an 



FIG. 2: Left panel: Nodal structure of superconducting gap. 
Here, we used the following parameters t = 4, a = 0, jx = 8, 
and At = 2. Middle panle: Winding number z^i, Eq. (40), as 



a function of surface momentum for the (001) face. The color 
scale is such that purple corresponds to z^i =0, whereas light 
gray corresponds to ui = +1. Right panel: Band structure 
on the (111) face. 



underlying crystal symmetry; it dictates the form of anti- 
symmetric spin-orbit coupling, and, subsequently, possi- 
ble pairing symmetries (at least within BCS theory). In 
the following, we will assume monoclinic crystal symme- 
try as an exampl e. ( For a possible relation of our model 
to BiPd, see Sec.[lV|. 

We start from the BdG Hamiltonian, 
H = {l/2)Ek[^{k)]^^{k)^{k), where ^{k) = 
{c^{k),c^{k),c\{-k),cl{-k)) and 



n{k) 



( m 

I At(fe) 



A{k) 
-h^{-k) 



(42) 



For NCSs under external magnetic field, we consider the 
following form of h{k) and A(/c), 

h{k) = e{k)l2x2 + al{k) - a^B-a, 
A{k) = (AJ2X2 + ApKk) ■ ^){icTy), (43) 

where e{k) is the band structure in the normal metal 
phase, a is the strength of spin-orbital coulping, l{k) rep- 
resents the orbital direction, a = (cri, (72, as) are Pauli 
matrices, B is the external (Zeeman) magnetic field, Ag 
and Ap are the gap function for singlet and triplet pair- 
ing. Here we assume the d-vector is locked in the same 
direction as the orbital direction l{k). For simplicity, we 
consider the normal state band structure in the cubic 
lattice, 

e{k) = t(cos kx + cos ky + cos kz) — /i, (44) 

where t is the hopping amplitude and /i is the chemical 
potential. The form of l{k) is set by the crystal symme- 
try. Here we consider C2 symmetry (monoclinic crystal 
symmetry), in which case l{k) can be written as ^ 



^ ai sin kx + sin ky ^ 

a2 sin ky + sin kx 
\ as sin kz J 



(45) 
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where a^^i^...^5 are a real parameter. Observe that we 
have the conditions, e{k) = e{—k) and l{k) = —l{—k). 
There are in total nine parameters /i, Ag, Ap and a^, 
z = 1,...,5. 

Particle-hole symmetry (P) for the BdG Hamiltonian 
with external magnetic field is 

(46) 



p-^U{k)P = -H*{-k), 



where P = ai I2x2- Time-reversal symmetry (T = 
UtIC, where JC is conjugate operator) for the BdG Hamil- 
tonian without external magnetic field {i.e. -B = 0) is 



U^^H{k)UT = 'H*{-k), 



(47) 



where Ut = ^2x2 ^ ^cr2- It is worth mentioning that as 
the external filed is tuning on, time-reversal symmetry is 
broken while particle-hole symmetry is still maintained. 

b. Nodal structures As an example, the parameters 
are chosen to be t = —1.0, /i = —2.0, a = 1, Ap = 1.0, 
henceforth, and we vary A5, ai = a2, as, = a^. The 
system symmetric is for exchanging x and y directions. 
Typical nodal lines that can be realized in this model are 
shown in in Fig.jSj In Fig.[3j one can observe all the nodal 
structures projected on t he (1 , 0, 0) and (0, 1, 0) surfaces 
are symmetric. For Figs. |3(c)| and |3(f)[ the nodal rings, 
when projected onto the (1, 0, 0) and (0,1,0) surfaces, 
do not overlap with each other, i.e. there will be the 
regions of non-zero topological charges that can lead to 
the surface flat band. 

In the following, we focus on th e t wo- an d four-nodal 
ring structures shown in Figs. 3(c)| and |3(f)[ respectively. 
For the two-ring case, the winding number is +2 (—2) 
for the loops encircling the nodal ring located at /c^ > 
{ky < 0) and is zero for the loops that do not encircling 
any nodal rings. For the four-ring case, the winding num- 
ber is +2 for the loops encircling the nodal rings located 
at {ky^kz) = (+, +) and (ky^k^) = (+,—), where the 
sign +/— indicates the value of momentum {ky and kz) 
is greater/less than zero. The winding number is —2 for 
the loops encircling the nodal rings at {ky^kz) = ( — ,+) 
and {ky^kz) = ( — ,—)• One can observe for this four-ring 
case, the winding numbers (topological charges) of each 
ring do not cancel out (when projecting on the (0,0,1) 
surface, say). Hence the flat bands can also appear at the 
(0, 0, 1) surface and will be discussed in the later section. 

c. Surface states The zero-energy flat bands can ap- 
pear in the surface Brillouin zone that come from the 
nodal rings projected on the surface. One can directly 
get these zero-energy flat bands from constructing the 
Hamiltonian in the position space with open boundary 
condition. In Fig. [4j the energy dispersion with respect 
to ky for open boundary conditio n in x direction is shown 
for the two nodal ring case (Fig. 3(b) ) at /c^ = 0, and for 
the four nodal ring case (Fig. 3(f7 ) at kz = 7r/2, respec- 
tively. One can see the flat dispersion for surface states 
appear for both two ring and four ring nodal structures. 
In both cases, their flat dispersion can contribute to the 
local density of states as we will show in the following sec- 
tion. The external magnetic field along z direction {i.e. 





(a) 



(b) 





(c) 



(d) 





(e) 



(f) 




FIG. 3: Typical nodal structures with different param- 
eters, (As,ai = a2,a3,a4 = as); (a)(1.0, 1.0, 0.5, 0.0) 
(b)(1.0, 0.9, 0.0, 0.1); (c) (1.0, 0.5, 0.0, 0.5) 

(d)(1.5, 0.7, 1.8, -0.7); (e)(1.0, 0.2, 1.1, -0.2) 

(f)(1.8, 0.6, 1.3, 0.6); (g)(0.6, 1.0, 0.0, 0.0). 



parallel to the surface) does not split o r shif t these fiat 
bands for the two nodal ring case (Fig. |4(b")] ). However, 
in the case where the external magnetic field along x di- 
rection {i.e. perpendicular to the surface), the flat bands 
get shifted to the opposite directions as shown in Fig. 
|4(c)[ (This is expected due to particle- hole symmetry). 

In Fig. [5| for the four nodal ring case, the flat band 
is located at (0,0, 1) surface within the nodal rings. As 
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the external magnetic filed turning on, the fiat band is 
shifted away from zero for the field both along z and x 
directions. 



The quasi-particle LDOS defined for BdG Hamiltonian is 




FIG. 4: Energy dispersion with open bounda ry in x direction 
corresponding to the nodal structure in Fig. 3(c)| at kz — 0, 
with 5 = (a), = 0.5 (b) and = 0.3 (c). (d) Energy 
dispersion with open boun dary i n x direction corresponding 
to the nodal structures Fig. 3(f) at A;^ = 7r/2 with B = 0. 




D{x,E) 



n,k± 



+ {\u^{x, k^)\^ + \u^{x, k^)\^)S{E + e(n, k^))] , (48) 

where x represents a site in x direction, 
{u^{k_\_)^ui{k±)^v^{kj_)^v^{k_\_)) is the n-th eigenvector 
with energy e{n,k±) for the Hamiltonian 1-L{x,x\k±) 
describing the Hamiltonian with fixed k±^ E is the 
energy and NyNz is the number of sites on the (1,0,0) 
surface. For the LDOS at the (1,0,0) surface, we set 
X = 1. 

In Fig. [6] for the case of two nodal rings, we plot the 
LDOS at surface in the presence of the external mag- 
netic field, B = (0,0, B^) and (5^,0,0). (We neglect 
the orbital effect, but focus on the effects of the Zeeman 
field.) As we increase 5, the zero energy peak is robust 
for the field along z direction but will be split for the 
field along x. Interestingly, for the former case, the zero 
energy peak grows up because that some of the states are 
pushed inward to the zero energy. 




i)3 0.0 0.5 
Energj 



FIG. 6: The quasiparticle LDOS at the (1, 0, 0) surface under 
external magnetic field for the case of two nodal rings, (a) 
B = (0,0,0), (b) (0,0,0.3), (c)(0, 0,0.5), and (d)(0.3,0,0) 
respectively. 



FIG. 5: Energy dispersion with open boundary at z with 
respect to four nodal rings structures in Fig. 3(f) at /c^ = 7r/2 



with (a) B = (0, 0, 0), (b) (0, 0, 0.1), (c) (0.3, 0, 0). 



The tunneling current is proportional to the local den- 
sity of states (LDOS) at surfaces. One should expect this 
fiat band will contribute to the LDOS to form a peak that 
is the origin of ZBCP for tunneling current measurement. 



D. Robustness of flat surface bands under disorder 

The topological protection of boundary modes appear- 
ing in nodal superconductors is weaker than the protec- 
tion of edge states in topological insulators or supercon- 
ductors with a full gap in the bulk. Nevertheless, as we 
explain below, the bound states in non-centrosymmetric 
superconductors enjoy a certain robustness, similar to 
the zero-energy edge states at the boundary of graphene 



11 



or dx'2-y2-vjme superconductors. Below we present three 
arguments suggesting that the surface states survive the 
presence of disorder. The first argument is that the topo- 
logical index can be defined even for a weakly disordered 
superconductor, namely by periodically repeating the fi- 
nite size system. The second and third arguments are 
based on semi-classical analysis and on symmetry con- 
siderations, respectively. 

(i) The bulk boundary correspondence for nodal super- 
conductors can be modified to incorporate dilute disor- 
der. This could be done, for example, by the well known 
method of periodically repeating the finite size system, 
which allows to define a winding number even in the pres- 
ence of disorder. [For recent theoretical work which has 
applied this method to study edge roughness in graphene 
see Refs. HIl [ZH] 

(ii) As has been discussed in Ref. 22, within a qua- 
siclassical description the zero-energy states in non- 
centrosymmetric superconductors occur whenever the or- 
der parameter on the negative helicity band changes sign 
along a quasiparticle trajectory. This criterion does not 
explicitly rely on the presence of translational symmetry, 
and hence suggests that the surface states are mostly un- 
affected by weak surface roughness or disorder (assuming 
that the semi-classical approximation still applies). 

(iii) As we have mentioned zero-energy boundary 
modes have been experimentally observed in graphene 
[73 - 75^ and in 6^3,2 -wave superconductors [76 - 78 , even 
though these systems are not necessarily clean. The infiu- 
ence of edge roughness and disorder on zero-energy edge 
states has been studied quite extensively in the literature, 
both for graphene [see, e.g., Refs. 79, 80 and d^'^-y^- 
wave superconductors [see, e.g., Refs. .81n83. ] The result 
of these investigations is that edge roughness and weak, 
dilute disorder does not remove the boundary states, but 
leads to a broadening of the energy spectrum and of the 
zero-bias conductance peak. 

To be more specific, it is known that, for example, 
in graphene the mixture of zigzag and armchair edges 
does not affect the existence of boundary modes at all. 
That is, in the limit of vanishing second neighbor hop- 
ping (t' = 0), the potential describing edge disorder in 
graphene does not have any matrix element that con- 
nects among the zero energy edge states [see Ref. [80] . 
This property is a consequence of the chiral symmetry 
(sublattice symmetry) of graphene. The source of insta- 
bility of the graphene edge states comes from the cou- 
pling between the bulk and the edge modes. Here how- 
ever, it is important to note that the bulk density of 
states near zero energy is small, and hence the coupling 
of the edge modes to the bulk is suppressed. For the 
case of the surface state of non-centrosymmetric super- 
conductors a similar argument can be repeated. Because 
of the presence of time-reversal and particle-hole symme- 
try, non-centrosymmetric superconductors also satisfy a 
chiral symmetry which protects the surface states against 
surface disorder. A detailed discussion of the infiuence of 
surface disorder on the zero-energy surface state of non- 



centrosymmetric superconductors will be the subject of 
a forthcoming publication. 

(iv) Finally, because of the chiral symmetry (sublattice 
symmetry) the surface disorder potential does not have 
any matrix element that connects among the zero energy 
edge states with positive (or negative) topological num- 
ber. For the study of the edge state in graphene, this 
argument has been mentioned in, e.g., Ref. [80] . 

The point is that nodal lines are only individually 
stable. But they are not robust against commensurate 
perturbations which hybridize nodal lines with opposite 
topological charge. Hence, one needs to assume the pres- 
ence of some sort of translational symmetry which pre- 
vents commensurate perturbations. However, this is not 
to say that all the translational symmetries of the lattice 
need to preserved. As explained above, edge roughness 
and weak dilute disorder in general does not lead to the 
hybridization of the nodal lines. Hence, this type of dis- 
order will not affect the surface bound states. 

To summarize, weak dilute disorder or surface rough- 
ness does not destroy the surface bound states. The bulk 
boundary correspondence discussed in the manuscript 
can be modified in a straightforward manner to include 
weak dilute disorder by periodically repeating the finite 
size system. It seems fairly obvious that weak dilute dis- 
order does not completely remove the surface states but 
just leads to a broadening of the energy spectrum. 



IV. DISCUSSION 

Based on previous works, we have developed a uni- 
fied, general description of topologically protected nodal 
point, line, and surfaces by taking into account the di- 
mensionalities as well as discrete symmetries. While 
fermion doubling and nesting should be properly taken 
into account, the local stability of fermi surfaces is gov- 
erned by a K-theoretic argument, in the way similar to 
the classification of bulk topological insulators and su- 
perconductors, and topological defects in insulators and 
superconductors. 

According to the bulk-boundary correspondence, one 
of the most direct signatures of the topological aspects 
of fermi surfaces are the existence of boundary modes. 
(In the case of superconductors with topologically pro- 
tected nodes, these are the Andreev bound states at 
an interface.) These can be probed by angle-resolved 
photoemission measurements, or by scanning tunnel- 
ing spectroscopy (STS) of the surface density of states 
(SDOS). Alternatively, reading the bulk-boundary dic- 
tionary backward, detection of the boundary modes al- 
lows one to access the information on the nodal structure 
in the bulk; This provides us a tool to map out supercon- 
ducting nodes, complementary to other methods such as 
angle-resolved thermal transport measurements. [84] 

We have illustrated these general principles in terms 
of simple examples in superconducting systems in sym- 
metry class AIII, and in class Dili in three dimen- 
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sions. While our theoretical model may be too simple 
(in that there is only one band per spin, say), one could 
learn general lessons on topological aspects of nodal non- 
centrosymmetric superconductors in general (see below). 

The topological properties of non-centrosymmetric su- 
perconductors (NCSs) have been studied theoretically 
and experimentally recently. j59l469^ Due to the absence 
of an inversion center in the crystal structure of NCSs, 
antisymmetric spin-orbit coupling leads to the admix- 
ture of spin-singlet and spin-triplet pairing states. When 
time-reversal symmetric, NCSs belong to symmetry class 
Dili, for which topologically stable nodal line can ex- 
ist. While there are various known NCSs, not all of 
them have nodal lines presumably due to small spin- 
orbital coupling, {e.g., ResW, MgioIrigBie, M03AI2C, 
and Re24Nb5). j62H65] Moreover, some of these NCSs 
also carry heavy fermion features {e.g., CePtsSi, Ulr, 
CeRhSia, and CelrSia) [66'-'69 which complicates the ex- 
perimental identification of the pairing symmetry due to 
the intricacies of strong electronic correlations. 

More recently, BiPd has been reported to be a su- 
perconductor below 3.87K. [59l |60] The spin-orbit cou- 
pling in BiPd is expected to be large. BiPd has a mon- 
oclinic crystal structure (P2i)[59, ^ at superconduct- 
ing states. From the directional point contact Andreev 
(PCAR) measurements, the zero-bias conductance peaks 
(ZBCP) have been observed for both crystallographic di- 
rections {I \\ b and I ±b with / is the current direction). 
Moreover, these ZBCPs for different surfaces do not split 
under the external magnetic field B which is set to be 
parallel to the current direction {B \\ I). 

While our model of the class Dili superconductor with 
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monoclinic crystal symmetry may be an oversimplifica- 
tion as a model of BiPd, it does share some properties 
which are well-compared with the above experiments. 
We demonstrate, with a certain choice of parameters, 
the fiat surface bands exist both for (1,0,0) and (0,0, 1) 
surfaces; This is due to relatively low symmetry of the 
crystal; For the case of higher crystal symmetry, e.g., 
if the system has cubic symmetry, as in the case of 
Li2Pda^Pt3_a;B, the fiat surface bands would be absent 
on both (1,0,0) and (0, 1,0) surfaces. These fiat surface 
bands contribute to a peak at zero bias. We also noted 
that this zero energy peak in the LDOS splits in general 
in the presence time-reversal symmetry breaking pertur- 
bations. Quite interestingly, the experiment in Ref. [60] 
does not show such splitting in the tunneling conductance 
under magnetic field. This may correspond to the situa- 
tions we found in Sec. |III C[ where time-reversal symme- 
try breaking does not immediately destroy the peak in 
the quasiparticle LDOS. 
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